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ABSTRACT  In  order  to  understand  GPS  time  comparison 
characteristics  associated  with  influences  from  SA  (Selective 
Availability)  effects,  a  number  of  experiments  were  carried  out  in 
early  1992.  The  data  sampling  time  associated  with  the  experiments 
was  10  seconds.  Satellite  tracking  time  was  1  hour.  The  results 
of  frequency  spectrum  analysis  clearly  show  that:  1)  Under  normal 
conditions,  GPS  time  comparison  spectrum  characteristics  are 
basically  level.  There  are  no  clear  periodic  phenomena.  Noise 
processes  are  phase  modulated  white  noise.  2)  When  GPS  time 
comparisions  are  influenced  by  SA  effects,  low  frequency  spectra 
are  relatively  clearer.  There  are  obvious  periodic  phenomena. 
Among  these,  5-15  minut  •  period  signals  account  for  approximately 
77%.  This  is  nothing  else  than  the  cause  which  leads  to  obvious 
drops  in  GPS  time  comparison  precision. 

KEY  WORDS  Global  positioning  system  Time  comparison 
Selective  availability  Spectrum  characteristic 
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1  INTRODUCTION 


Since  November  15th  1991,  the  majority  of  BLOCK  II  satellites 
have  turned  from  normal  conditions  to  SA  configuration  operations. 
BLOCK  I  satellites  operate  under  normal  conditions  from  beginning 
to  end.  SA  (selective  availability)  conditions  are  a  number  of 
parameters  added  into  random  modulation  with  regard  to  signal 
characteristics.  How  about  GPS  time  comparisons  influenced  by  SA 
effects  with  respect  to  characteristics?  This  is  a  particular 
concern  for  GPS  time  transmission  users.  The  reason  is  that  the 
accuracy  and  precision  of  GPS  time  comparisions  influenced  by  SA 
both  drop  to  different  degrees.  There  are  some  results  which  fall 
relatively  greatly.  This  will  severely  influence  long  range  GPS 
time  synchronicity  results.  For  this  reason,  comparisons  and 
analyses  are  carried  out  from  received  data  with  regard  to  the 
results  associated  with  the  two  conditions.  In  conjunction  with 
this,  thinking  up  ways  to  reduce  the  influences  of  SA  effects  has 
already  turned  into  an  important  topic  associated  with  GPS  time 
transmissions . 

For  the  sake  of  comparing  GPS  time  comparison  results  under 
the  two  types  of  conditions,  we  carried  out  a  series  of  test 
measurements  and  analyses  in  early  1992.  Variances  associated  with 
time  domain  characteristic  analyses  clearly  showed  that  GPS  time 
comparison  accuracies  influenced  by  SA  effects  (use  was  made  of 
recieved  data  for  which  the  sampling  period  was  10  seconds  within 
a  tracking  period  of  13  minutes)  dropped  5-10  fold  as  compared  to 
the  results  for  normal  conditions.  This  article  primarily 
describes  analysis  results  for  a  number  of  frequency  domains 
associated  with  GPS  time  comparisions  under  the  two  types  of 
conditions.  From  these,  it  is  possible  to  know  a  number  of 
characteristics  and  patterns  of  SA  effects  influencing  GPS  time 
comparison  fequency  domain  characteristics. 

2  FREQUENCY  DOMAIN  CHARACTERISTIC  ANALYSIS 

Fig.l  is  an  example  of  changes  in  time  difference  measurement 
values  associated  with  GPS  time  comparisons  under  the  two 
conditions.  In  Fig.l,  comparison  signals  are  [UTC(SO) -GPS (PRN  J) ] 
(J  is  the  satellite  PRN  serial  no.).  Data  sampling  periods  are  10 
seconds.  Satellite  tracking  periods  are  1  hour.  The  situation 
associated  with  BLOCK  I  satellites  is  similar  to  the  situation 
associated  with  PRN  3  in  Fig.'s.  Amplitude  fluctuations  are  not 
great.  The  case  associated  with  BOCK  II  satellites  influenced  by 
SA  effects  is  similar  to  the  situation  associated  with  PRN  16  in 
Fig.'s.  Amplitude  fluctuations  appear  as  clear  wave  forms. 

The  purpose  of  frequency  domain  characteristic  analyses  lies 
in  precisely  specifying  whether  or  not  obvious  low  frequency  noise 
or  periodic  phenomena  exist.  It  is  possible — on  the  foundation  of 
time  domain  autocorrelation  function  analyses — to  go  through 
characteristic  analyses  associated  with  power  spectrum  density 
functions  and  periodograms  in  order  to  obtain  them. 
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Power  spectrum  density  (PSD)  functions  are  Fourier  transforms 
of  autocorrelation  functions.  Making  use  of  autocorrelation 
function  Rx(k)  obtained  from  GPS  time  comparison  dispersion  data 
X(t),  it  is  possible  to  calculate  power  spectrum  density  Sx(f)  to 
be: 


N-  l 

S*(/)  =  2[Ri(0)  +2  Rx(k) cosink f]t  f  =  j/Ti  0'  =  1,  2,  •••,  N/ 2) 

i-i 


(1) 


In  this.  N  is  the  number  of  data  sampling  iterations.  T  is  the 
time  interval  associated  with  finite  data.  Sx(f)  can  be  calcula  e 
by  FFT  algorithm  directly  from  GPS  time  comparison  equations  of 
time  X(t) .  Sx(f)  calculation  results  make  use  of  spectrum  windows 
(for  example,  Hamming  or  Hann  spectrum  windows)  to  carry  out 
smoothing  in  order  to  reduce  the  level  of  dispersion.  /50 


Fig.l  Time  Difference  Changes  Associated  with  GPS  Time 
Comparisons  Under  the  Two  Conditions 
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Fig. 2  Correlation  Diagrams  and  Frequency  Spectrum  Diagrams 
Associated  with  GPS  Time  Comparisons  Under  the  Two  Conditions 
Key:  (1)  Correlation  Diagrams  (2)  Frequency  Spectrum  Diagrams 


From  the  state  of  changes  in  power  spectrum  densities  Sx(f)  as 
a  function  of  Fourier  frequencies  f,  it  is  possible  to  precisely 
determine  a  number  of  frequency  domain  characteristics  associated 
with  variables  X(t):  (1)  If  this  type  of  change  is  basically 

level  or  undulates  very  little,  X(t)  is  a  modulated  phase  white 
noise  process.  (2)  If  the  slope  associated  with  this  type  of 
change  presents  an  exponential  attenuation,  X(t)  is  a  low  frequency 
noise  process.  (3)  If  this  type  of  change  has  clear  peak  values, 
X (t)  manifests  obvious  periodic  phenomena.  If  one  needs  to 
precisely  specify  periodic  phenomena,  it  is  possible  to  go  through 
Fourier  transforms  to  carry  out  periodogram  analyses. 

The  Fourier  expression  associated  with  dispersion  variable 
X(t)  is: 


X(0  =  «0  + 


X]  (  GljCOS 

/-l 


2«  jt 

N 


bjSin 


2  sr/t 

N 


) 


(2) 
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Table  1  Analysis  Results  for  BLOCK  I  Satellite  GPS  Time  Comparison 
Frequency  Domain  Characteristics  and  Time  Domain 
Characteristics 
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Key:  (1)  Reception  Satellite  (2)  Reception  Date  and  Start  Time 

(3)  Spectrum  Density  Function  Characteristics  (4)  Period 
Characteristics  (5)  Autocorrelation  Function  Characteristics 
(6)  Noise  Process  Slopes  (7)  Basically  Level  (8)  Not  Clear 
(9)  Basically  Level 
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In  equations, 

k=\N/- 
kN-  l)/2 

t  =  1,2,—,  Nt  Aj  =  v'  0/  ±bY' 


uh 


0  j  =  arcsin(a,/  A ,) 


By  FFT  methods,  solve  for  parameters  ao,  a j ,  and  b j .  After  that, 
incorporate  them  into  amplitude  Aj  and -phase  .  Carrying  out  F 
tests  on  the  basis  of  these  parameters,  it  is  possible  to  precisely 
determine  whether  or  not  certain  periodic  terms  are  obvious 
(P=l/f=T/j) .  During  inspections,  obvious  levels  are  selected  as 
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a=0 . 01.  In  conjunction  with  this,  the  previous  3  obvious  periods 
are  selected  to  do  comparisons.  /51 

Table  2  Analysis  Results  for  BLOCK  II  Satellite  GPS  Time 
Comparison  Frequency  Domain  Characteristics  and  Time  Domain 
Characteristics 
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Key:  (1)  Reception  Satellite  (2)  Reception  Date  and 

Start  Time  (3)  Spectrum  Density  Function  Characteristics  (4) 
Period  Characteristics  (Arranged  in  accordance  with  the  order  of  ■ 
obviousness.  Time  unit:  min)  (5)  Autocorrelation  Function 
Characteristics  (6)  Noise  Process  Slope  (7)  Low  Frequency 
Spectrum  Relatively  Clear  (8)  Basically  Level  (9)  Wave  Shaped 
Variation 


6 


Results  associated  with  frequency  domain  characteristics  and 
time  domain  characteristics  are  corresponding  ones.  In  the  Fig.'s 
and  Tables,  we  take  their  results  and  set  them  out  together  to  make 
comparisons.  Fig. 2  is  examples  of  autocorrelation  diagrams  and 
frequency  spectrum  diagrams  associated  with  GPS  time  comparisons 
under  the  two  types  of  conditions.  Table  1  and  Table  2  are, 
respectively,  analysis  results  for  frequency  domain  characteristics 
and  time  domain  charcteristics  associated  with  GPS  time  comparisons 
for  BLOCK  I  satellites  and  BLOCK  II  satellites.  In  these  two 
tables,  time  domain  results  have  autocorrelation  function 
characteristics  and  noise  process  slopes  (drawn  up  in  accordance 
with  different  sampling  period  ranges  and  two  types  of  Allan 
variances,  the  majority  of  autocorrelation  coefficients  are  larger 
than  0.9) [1,2].  As  far  as  the  explanations  of  these  Fig.'s  and 
tables  are  concerned,  when  sampling  periods  are  1  minute  and 
satellite  tracking  periods  are  1  hour,  1)  frequency  spectrum 
characteristics  associated  with  GPS  time  comparisons  for  normal 
conditions  (all  BLOCK  I  satellites  and  BLOCK  II  satellites)  are 
basically  level.  Periodic  characteristics  are  not  clear  and  are 
phase  modulated  white  noise  processes.  2)  Low  frequency  spectra 
associated  with  GPS  time  comparisons  influenced  by  SA  effects  (a 
majority  of  BLOCK  II  satellites)  are  relatively  clear.  Periodic 
phenomena  are  obvious.  In  the  preceeding  3  obvious  periods,  5-15 
minute  period  signals  account  for  77%.  This  is  nothing  else  than 
the  cause  for  drops  in  precision  associated  with  GPS  time 
comparisons  influenced  by  SA  effects.  This  also  explains  why  there 
are  some  statistical  patterns  associated  with  it. 
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